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We carried out observations of the central and 20' east offset regions of the cluster of galaxies Abell 1060 
with Suzaku. Spatially resolved X-ray spectral analysis has revealed temperature and abundance profiles 
of Abell 1060 out to 27' ~ 380 hj^ kpc, which corresponded to ~ 0.25 ttso- Temperature decrease of the 
intra cluster medium from 3.4 keV at the center to 2.2 keV in the outskirt region are clearly observed. 
■ Abundances of Si, S and Fe also decrease by more than 50% from the center to the outer, while Mg shows 

fairly constant abundance distribution at ~ 0.7 solar within r 17'. O shows lower abundance of ~ 0.3 
solar in the central region (r < 6'), and indicates a similar feature with Mg, however it is sensitive to 
the estimated contribution of the Galactic components of kT\ ~ 0.15 keV and kTi ~ 0.7 keV in the outer 
annuli (r > 13'). Systematic effects due to the point spread function tails, contamination on the XIS filters, 
instrumental background, cosmic and/or Galactic X-ray background, and the assumed solar abundance 
tables are carefully examined. Results on temperature and abundances of Si, S, and Fe are consistent with 
those derived by XMM-Newton at r < 13'. Formation and metal enrichment process of the cluster are 
discussed based on the present results. 
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1. Introduction 

Clusters of galaxies, being the largest virialized system 
in the universe, are filled with the intracluster medium 
(ICM), which consists of X-ray emitting hot plasma with 
typical temperatures of a few times 10 7 K. X-ray spec- 
troscopy of the ICM can immediately determine its tem- 
perature and metal abundances. The metal abundances 
of ICM have a lot of information to understand the chem- 
ical history and evolution of clusters. A large amount 
of metals of the ICM are mainly produced by supernovae 
(SN) in early-type galaxies (Arnaud et al. 1992; Renzini et 
al. 1993), which are classified roughly as Type la (SN la) 
and Type II (SN II). Si, S and Fe are synthesized in both 
SN la and SN II, while a elements such as O, Ne, and Mg 
are mainly in SN II, which are explosions of massive stars 



with initial mass above ~ 10 Mq. The metals produced in 
the galaxies are transfered into the ICM by galactic wind 
and/or ram pressure stripping. 

ASCA firstly revealed the distribution of Si and Fe 
in the ICM (Fukazawa et al. 1998; Fukazawa et al. 
2000; Finoguenov et al. 2000; Finoguenov et al. 2001). 
The derived iron-mass-to-light ratios (IMLR; Ciotti et 
al. 1991; Renzini et al. 1993; Renzini 1997) are nearly 
constant in rich clusters and decrease toward poorer sys- 
tems (Makishima et al. 2001). Recent observations with 
Chandra and XMM-Newton allowed detailed studies of 
the metals in the ICM. These observations, however, 
showed abundance profiles of O, Mg, Si and Fe only for 
the central regions of very bright clusters or groups of 
galaxies dominated by cD galaxies in a reliable manner 
(Finoguenov et al. 2002; Fukazawa et al. 2004; Matsushita 
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ct al. 2003; Tamura et al. 2003). The abundance pro- 
files of O and Mg, in particular for the cluster outer re- 
gions, are still poorly determined, because these satel- 
lites are characterized by relatively high intrinsic back- 
ground levels. Tamura et al. (2004) derived IMLR for 
five clusters within 250 h^ kpc to be <~ 0.01 M Q /i Q , 
and the oxygen mass within 50 kpc for several clus- 
ters. However, oxygen-mass-to-light ratios (OMLR) for 
rich clusters are not reliable due to the lower emissivity of 
O VII and O VIII lines in higher temperatures. De Grandi 
& Molendi (2001); Hayakawaet al. (2006) found that clus- 
ters associated with cD galaxies and central cool compo- 
nents showed abundance concentration in the cluster cen- 
ter, while clusters without cD galaxies showed flatter pro- 
files. The central metallicity enhancement in the cool core 
clusters were further studied and the excess metals were 
shown to be supplied from the cD galaxies (De Grandi et 
al. 2004). 

Abell 1060 (hereafter A 1060) is a nearby cluster of 
galaxies (z — 0.0114) characterized by a smooth and sym- 
metric distribution of intraclustcr medium (ICM) , and has 
no cD-galaxy at the center. Chandra observation detected 
very compact X-ray emissions from the central two ellip- 
tical galaxies (Yamasaki et al. 2002), and the tempera- 
ture and abundance distributions were shown to be some- 
what inhomogeneous at the cluster center (Hayakawa et 
al. 2004). Hayakawa et al. (2006) detected a temperature 
drop by ~ 30% from the central region to r ~ 13' with 
XMM-Newton, while it had previously been considered as 
flat on the basis of the ASCA and ROSAT observations 
(Tamura ct al. 1996, Furusho et al. 2001). 

This paper reports results from Suzaku observations 
of A 1060. Owing to the low-background nature of the 
Suzaku XIS, we are able to measure the temperature and 
abundance profiles to a much outer region than the previ- 
ous XMM-Newton study. We use H = 70 km s" 1 Mpc -1 , 
Cl\ = 1 — S1m = 0.73 in this paper. At a redshift of 
z = 0.0114, 1' corresp onds to 14 kp c, and the virial ra- 
dius, r-180 = 1-95 h^ o ^k{T}/l0 keV Mpc (Markevitch ct 
al. 1998), is 1.53 Mpc for the average temperature of 
(T) = 3 keV (Hayakawa et al. 2006). Throughout this 
paper we adopt the Galactic hydrogen column density of 
A H = 4.9 x 10 20 cm" 2 (Dickey & Lockman 1990) in the di- 
rection of A 1060. Otherwise noted, the solar abundance 
table is given by Anders & Grevesse (1989), and errors are 
90% confidence region for a single interesting parameter. 

2. Suzaku Observation and Data Reduction 

2.1. Observation 

Suzaku carried out two pointing observations for A 1060 
in November 2005, the central region and 20' east offset 
region, with exposure of 40.5 and 55.5 ks, respectively. 
The observation log is summarized in table 1, and the 
combined X-ray Imaging Spectrometers (XIS; Koyama 
ct al. 2006) image in the 0.5-7 keV range is shown in 
figure 1. We utilize only the XIS data in this paper. 
The XIS is an X-ray CCD camera, which consists of 
one back- illuminated sensor (BI = XIS1) and three front- 



Table 2. Estimated column density of the contaminant for 
each sensor at the center of CCD in unit of 10 18 cm~ 2 . 





XIS0 


XIS1 


XIS2 


XIS3 


Carbon 


1.37 


1.83 


2.52 


4.04 


Oxygen 


. ... 0.228 


0.305 


0.420 


0.673 



illuminated sensors (FI = XIS0, XIS2, XIS3). The BI 
and FI sensors have different advantages. The former 
has higher quantum efficiency in the soft energy band 
(2? < 1 keV), while the latter shows lower instrumental 
non X-ray background (NXB). The XIS was operated in 
the Normal clocking mode (no window nor burst option, 
so 8 s exposure per frame) , with the standard 5 x 5 or 3 x 3 
editing mode (Koyama et al. 2006). 

The optical blocking filters (OBF) of the XIS have been 
gradually contaminated in time by out-gas from the satel- 
lite, and the degradation of the low energy transmission 
was already significant in November 2005. The thickness 
of the contaminant is different among sensors, and is also 
dependent on the location on the CCD. The estimated 
column density (C/0=6 in number ratio is assumed) at 
the center of the CCD is listed in table 2, and the cal- 
culated X-ray transmission for BI (XIS1) at each annular 
region is plotted in figure 2. 1 Though the thickness of the 
OBF contaminant is different on each annulus, this effect 
is considered in the calculation of the Ancillary Response 
File (ARF) by the "xissimarfgen" Ftools task (Ishisaki et 
al. 2007). The energy resolution was also degraded slightly 
(FWHM - 150 eV at 5.9 keV) after the launch, due to the 
radiation damage of the CCD. 

2.2. Data Reduction 

We used the version 0.7 processing data (Mitsuda et 
al. 2007), and the analysis was performed with HEAsoft 
version 6.0.6 and XSPEC 11.3.2t. We started the event 
screening from the cleaned event file, in which selection 
of the event grade and bad CCD column, disposal of non- 
observational intervals (during maneuver, data-rate low, 
South Atlantic Anomaly, Earth occultation, bright Earth 
rim to avoid scattered solar X-ray), and removal of hot 
and flickering pixels by the "cleansis" Ftools, were already 
conducted. The exposure time given in table 1 is for the 
cleaned event file. We further applied the Good-Time 
Intervals (GTI) given for excluding the telemetry satura- 
tion by the XIS team. The light curve of each sensor in 
the 0.3-10 keV range with 16 s time bin was also examined 
to reject periods of anomalous event rate greater or less 
than ±3cr around the mean. After the above screenings, 
remaining exposure of the central observation was 40.2 ks 
and that of the offset observation was 52.9 or 54.6 ks (FI 
or BI) . These exposures are not so different from those in 
table 1, which represent that the NXB was almost stable 
during the both observations. The event screening with 
the cut-off rigidity (COR) was not performed in our data. 

1 The calibration database file of ae_xiW_contami_20060525.f its 

was used for the estimation of the XIS contamination (N = 
0,1,2,3 corresponding to the XIS sensor). 
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Fig. 1. (a) Combined XIS image of the central and offset observations in the 0.5-7.0 keV energy range. The observed XISO-3 images 
are added on the sky coordinate after removing each calibration source region, and smoothed with a = 16 pixel ~ 17" Gaussian. 
Estimated components of cxtragalactic X-ray background (CXB) and instrumental background (NXB) are subtracted, and exposure 
and vignetting are corrected, (b) Exposure map in ks, in which coordinates of the image are the same with (a). Four XIS sensors 
are treated separately. Small exposure at the CCD corners correspond to the 55 Fe calibration source locations. The annular regions 
utilized in the spectral analysis are indicated by green circles. 



Table 1. Suzaku observation of A 1060 and NGC 2992 



Target name Sequence number Date Exposure time (RA, Dec) in J2000 * (I, b) 



A 1060 center 
A 1060 offset 


800003010 
800004010 


2005-Nov-22 
2005-Nov-20 


40.5 ks 
55.5 ks 


(10 h 36 m 42?8, 
(10 h 38 m 03?8, 


-27°31'42") 
-27°31'42") 


(269.°60, 26.°49) 
(269.°88, 26P65) 


NGC 2992 T 
NGC 2992 t 
NGC 2992 1 


700005010 
700005020 
700005030 


2005-Nov-06 
2005-Nov-19 
2005-Dec-13 


38.8 ks 
38.7 ks 
48.3 ks 


(09 h 45 m 45?2, 
(09 h 45 m 41?7, 
(09 h 45 m 51f0, 


-14°16'10") 
-14°16'05") 
-14°16'51") 


(249.°69, 28.°83) 
(249.°66, 28.°82) 
(249.° 70, 28.° 84) 



* Average pointing direction of the XIS, written in the RA_NOM and DECJMOM keywords of the event FITS files, 
t We used NGC 2992 data to estimate the Galactic component in subsection 4.5. 
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Fig. 2. (a) Estimated transmission of the contaminant on the BI (XIS1) sensor for each annular region used in the spectral analysis 
plotted against the X-ray energy. Transmissions in the central observations are drawn by solid lines, and the offset observations are 
by dashed lines. These transmissions are calculated by "xissimarfgen" with a calibration file of ae_xil_contami_20060525 . f its, and 
written in the CONTAMI.TRANSMIS column of the ARF response, (b) The calculated transmission plotted against each annular region 
in the energies of 0.5, 0.6, 0.7, and 1.0 keV. 
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2.3. NXB & CXB Subtraction 

In order to subtract the NXB and the extra-galactic cos- 
mic X-ray background (CXB; see Brandt & Hasinger 2005 
for review), we used the night earth database of 770 ks 
exposure provided by the XIS team for the NXB, and es- 
timated the CXB component using the ASCA results. 

The night-earth spectra were extracted from the same 
detector region as the A 1060 observation in order to can- 
cel positional variation of the NXB. Furthermore, we di- 
vided the night-earth data referring to COR in the ranges 
of < 4 GV, 4-13 GV in 1 GV step, and > 13 GV, be- 
cause it is known that the intensity and energy spectrum 
of the NXB are primarily correlated with COR at the 
orbital location of the satellite. The NXB spectra to be 
subtracted in the spectral fitting were estimated by adding 
these COR-sorted night-earth spectra weighted with expo- 
sure times of the A 1060 observation in the corresponding 
COR range. 

The CXB component was estimated by the "fake" com- 
mand of XSPEC using uniform-sky ARFs, which is gener- 
ated by "xissimarfgen" assuming that the whole sky (prac- 
tically, r < 20') has the constant intensity and the same X- 
ray spectrum. We assumed a power-law spectrum for the 
CXB with the values by Kushino et al. (2002), r = 1.4 and 
Sx = 5.97 x 10" 8 erg cm" 2 s" 1 sr" 1 (2-10 keV), absorbed 
with the neutral hydrogen column of Ah = 4.9 x 10 20 cm~ 2 
(Dickey & Lockman 1990). The above CXB intensity is 
taken from table 3 of Kushino et al. (2002), for the in- 
tegrated spectrum with source elimination brighter than 
S a = 2 x 10~ 13 erg cm- 2 s" 1 (2-10 keV) in the GIS filed of 
view with r = 1.4 (fix) and the nominal NXB level (0%). 
It is confirmed that this gives a reasonable estimate of the 
CXB contribution for the XIS in subsection 6.2 of Ishisaki 
et al. (2007) and Fujimoto et al. (2007) 

Figures 3 (a)-(c) show the background level for Suzaku 
BI, FI sensors, and XMM-Newton MOS1, at the same 6- 
9' annulus of the cluster. The observed spectrum after 
the CXB and NXB subtraction are compared with the 
estimated CXB and NXB spectra separately for Suzaku, 
while sum of the CXB and NXB are indicated for the 
XMM-Newton because the blank-sky data was used as the 
background. These figures show the background level of 
Suzaku XIS is ~ 50% lower than that of XMM-Newton. 
It is also notable that there is a strong A1-K Q peak at 
1.49 keV for the XMM-Newton background, which makes 
the determination of the Mg abundance quite difficult. 
The S/N ratio of the Suzaku XIS BI/FI sensors are by 
about 1.8/2.0 times higher than that of the XMM-Newton 
MOS at 1 keV, and 1.3/1.8 times higher at 4 keV. The 
particle background is more stable for Suzaku than XMM- 
Newton due to its low-earth orbit. 

3. XMM-Newton Observation &: Analysis 

The XMM-Newton observation of A 1060 was carried 
out on 2004 June 29 with total exposure of 64 ks. The 
same dataset with Hayakawa et al. (2006) was utilized. We 
used only the MOS data, and data reduction and analy- 




Projected radius (') 



Fig. 4. In the upper panel, a radial profile of the sur- 
face brightness of A 1060 in the 0.8-3 keV band are plot- 
ted for XMM-Newton MOS1+2 (r < 13'). The best-fit 
double-/3 model is shown by the solid gray line, and the two 
/3-components are indicated by dashed lines. In the bottom 
panel, fit residuals are shown in unit of a. 



Table 3. The best-fit parameters of figure 4. 








r c 


q * 

^0.8-3 keV 


narrower . . . 


. 1.3 ±0.1 


3.7' ±0.3' 


1.26 ±0.14 


wider t .... 


. 0.69 (fix) 


7.3' (fix) 


0.96 ±0.05 



* "So. g-3 kov at the center in unit of 10 7 counts s 1 arcsec 2 . 
' /3 and r c is fixed to the values by Hayakawa et al. (2006). 



sis were performed with SAS version 6.0 and HEAsoft 
version 6.0.6. As for the XMM-Newton data, data re- 
duction and spectral analysis were based on Sato et al. 

(2005) and Hayakawa et al. (2006). Hayakawa et al. (2006) 
have reported the temperature and abundance profiles 
with XMM-Newton observation, however, we reanalyzed 
the data in order to compare with the Suzaku results on 
the same criterion. The differences form Hayakawa et al. 

(2006) are: (1) the event extraction region, (2) using the 
vapec thin thermal plasma model (Smith et al. 2001) in- 
stead of vmekal (Mewe et al. 1985; Mewe et al. 1986), and 
(3) deriving the element abundances separately. 

Utilizing the XMM-Newton image, we also derived the 
surface brightness profile of A 1060 needed to generate the 
Suzaku ARFs, because the spatial resolution of XMM- 
Newton is superior to that of Suzaku. Figure 4 shows 
a radial profile of A 1060 in the 0.8-3 keV energy range, 



Temperature and Abundance Profiles of Abell 1060 with Suzaku 



■5 




(<0 



8'<r<9' 1 




Energy (keV) 



Energy (keV) 



1 2 
Energy (keV) 



Fig. 3. (a) The observed spectrum at the annular region of 6-9' with the Suzaku BI (XIS1) sensor is plotted in black crosses after 
subtracting the estimated CXB and NXB components, which are plotted by gray crosses and a black histogram, respectively, (b) 
Same as (a) but for the FI (XIS0+XIS2+XIS3) sensors. The count rate drop below < 0.4 keV is due to the event threshold of the FI 
CCDs (Koyama et al. 2006). Although 0.2-10 keV energy range is shown here, only 0.4-7.1 keV band was utilized for the spectral 
fitting for both BI and FI. (c) Same as (a) but for the XMM-Newton MOS1 sensor, and the estimated CXB+NXB component using 
blank-sky observations are plotted in gray. The 0.5-8.0 keV energy range was used for the spectral fitting with XMM-Newton. 



fitted with a double-/? model. The origin of the radial pro- 
file is placed at (RA, Dec) = (10 h 36 m 43fl, -27°31'46") in 
J2000. The best-fit parameters are summarized in table 3, 
in which (3 and r c of the wider /3-component are fixed to 
the values by Hayakawa et al. (2006). 



Table 4. Area, coverage of whole annulus, S0URCE_RATI0_ 
REG and observed counts for each annular region. SOURCE. 
RATlO_REG represents the flux ratio in the assumed spatial 
distribution on the sky (double-/? model) inside the accumu- 
lation region to the entire model, and written in the header 
keyword of the calculated ARF response by "xissimarfgen" . 



4. Spectral Analysis and Results 

4-1. Suzaku XIS Spectra 

We extracted spectra from seven annular regions of 0- 
2', 2-4', 4-6', 6-9', 9-13', 13-17' and 17-27', centered 
on (RA, Dec) = (10 h 36 m 42f8, -27°31'42"). The first 
four annuli were taken from the central observation, and 
the rest of thee annuli were from the offset observation. 
Table 4 lists areas of the extraction regions (arcmin 2 ), cov- 
erage of the whole annulus (%), the SOURCE_ratio_REG 
values (%; see caption for its definition) and the observed 
counts in 0.4-7.1 keV including NXB and CXB for the BI 
and FI sensors. 

Although the thickness of the OBF contamination is dif- 
ferent among sensors as shown in table 2 and figure 5(b), 
we have confirmed that the four sensors give quite con- 
sistent fit results after incorporating the contamination 
effect into the ARFs, as demonstrated in figure 5(a). We 
therefore add the three FI spectra (XIS0, XIS2, XIS3), 
hereafter. 

Each annular spectrum is shown in figure 6. The ionized 
Mg, Si, S, Fe lines are clearly seen in each ring. The O VII 
and O VIII lines are prominent in the outer rings, however, 
most of the OVII emission is supposed to come from the 
local Galactic emission, which will be examined in detail 
in subsections 4.3-4.5. 

4-2. Strategy of Spectral Fit 

The basic strategy of the spectral fit is described in sub- 
section 6.3 of Ishisaki et al. (2007). The observed spec- 
trum is assumed to contain (a) thin thermal plasma emis- 
sion from the ICM, (b) local Galactic emission, (c) CXB, 
and (d) NXB. The estimation of (c) and (d) is described 
in subsection 2.3, and both are subtracted from the ob- 
served spectrum. The spectrum of (b) can be represented 
by one or two thin-thermal plasma model(s), apec, with 



Region 


* Area T 


Coverage 


f SOURCE. 


* Counts 5 




(arcmin 2 




RATIO_REG 


BI 


FI 


0-2' 


12.6 


100.0% 


7.7% 


30,549 


64,606 


2-4' 


37.7 


100.0% 


13.7% 


55,186 


114,220 


4-6' 


62.8 


100.0% 


13.8% 


49,259 


100,080 


6-9' 


137.0 


96.9% 


17.2% 


52,719 


100,750 


9-13' 


49.1 


17.8% 


2.8% 


10,696 


17,329 


13-17' 


76.2 


20.2% 


2.2% 


11,032 


21,010 


17-27' 


167.3 


12.1% 


2.0% 


18,055 


28,834 



* The first four annuli arc extracted from the central observa- 
tion, and others are from the offset observation. 

t The largest values among four sensors are presented. 

* sOURCE_RATlO_REG = Coverage J r r ° ut S(r)rdr/ f^ 1 "^ S(r)rdr, 

where S(r) represents the radial profile shown in figure 4 and 
table 3, and the r max is set to 30'. 

§ Observed counts including NXB and CXB in 0.4-7.1 keV. 

1 solar abundance, however, it may vary from field to field 
by more than an order of magnitude (Kushino et al. 2002). 
The ICM spectrum (a) can be represented by a variable 
abundance thin-thermal plasma model, vapec, whose best- 
fit parameters are what we want. An important point is 
that the spatial distributions are different between (a) and 
(b). The former follows the surface brightness of the clus- 
ter, while the latter is supposed to have almost a uniform 
distribution in the XIS field of view. 

We therefore generated two different ARFs for the spec- 
trum of each annulus, A v and A B , which respectively 
assume the uniform-sky emission and ~ 1° x 1° size of 
the double-/? surface brightness profile obtained with the 
XMM-Newton data (table 3). A u was used to evaluate the 
CXB, and the surface brightness of the Galactic compo- 
nent (b) at each annulus in combination with the XSPEC 
"fakeit" command, and A B was used for the actual fitting. 
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Fig. 5. (a) The uppermost panel show the observed spectra after subtracting the estimated CXB and NXB components at the 
annular region of 2—4' plotted separately for the XISO (black), XIS1 (red), XIS2 (green), and XIS3 (blue) sensors. The cross markers 
denote the observed spectra and the solid lines show the best-fit model with apec 1 + apec 2 + phabs X vapec. The energy range around 
the Si K-edge (1.825-1.840 keV) is ignored for the spectral fit. Lower three panels show the residuals of the fit in unit of a for XIS1, 
XISO, XIS2, and XIS3 from upper to lower, respectively, (b) Plots of the calculated XIS effective area (ARF + RMF) including the 
XIS quantum efficiency. These responses are used in the spectral fit in (a) with the same colors. The quantum efficiency is much 
higher for BI (XIS1; red) than FI sensors, while the OBF contamination is thicker in the order of XIS3, XIS2, XIS1, and XISO. 
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Fig. 6. The upper panels of (a)-(g) show the background (NXB+CXB) subtracted 
spectra at the annular regions which is denoted in the panels, and they are plot- 
ted by red and black crosses for BI and FI, respectively. The estimated CXB 
and NXB components are subtracted, and the plotted data are fitted with the 
apec 1 + apec 2 + phabs X vapec model drawn by green and yellow lines for the BI and 
FI spectra. The apec 1 and apec 2 components representing the local Galactic emission 
for the BI spectra are indicated by cyan and orange lines. The energy range around 
the Si K-edge (1.825-1.840 keV) is ignored for the spectral fit. The energies of the 
several prominent lines are also indicated in the panels. The lower panels show the 
fit residuals in unit of a. 
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We have confirmed that the assumed double-/? surface 
brightness profile is consistent with the observed Suzaku 
image by about ±15% in figure 5 of fshisaki et al. (2007). 
We further constrained that the surface brightness of the 
Galactic component (b) is nearly constant among all the 
extraction annuli, whereas its spectral shape is determined 
from the spectral fit of our Suzaku data. Details will be 
described in the next subsection. 

We adopted the nominal Redistribution Matrix Files 
(RMF) of ae_xiA_20060213.rmf for spectral fitting, al- 
though slight degradation in energy resolution is expected 
(subsection 2.1). The ARFs were generated by "xissi- 
marfgen" , and were convolved with the RMFs and added 
for three FI sensors, using the "marfrmf" and "addrmf" 
tasks in Ftools. The spectra from BI and FI are fitted 
simultaneously in the 0.4-7.1 keV band except for energy 
range of anomalous response around the Si K-edge (1.825- 
1.840 keV). We ignored below 0.4 keV because the C edge 
(0.284 keV) seen in the BI spectra could not be reproduced 
perfectly in our data. Energy range above 7.1 keV was also 
ignored because background Ni line (~ 7.5 keV) left arti- 
ficial structures after the NXB subtraction at large radii. 
It is also known that the XIS response in E > 8 keV are 
not fully understood at the present stage. In the simulta- 
neous fit of BI and FI, only the normalization are allowed 
to be different between them, although we found that the 
derived normalizations are quite consistent between the 
two. 

4--S. Estimation of Galactic Component 

We found that the estimation of the Galactic compo- 
nent significantly affect the determination of the oxygen 
abundance, because both the ICM and the Galactic com- 
ponent contribute to the O VIII emission line (0.653 keV), 
and the XIS cannot resolve them by redshift due to the 
limited energy resolution. On the other hand, most of the 
OVII emission is supposed to originate in the Galactic 
component, because the ICM temperature kT > 2 keV 
~ 2.3 x 10 7 K is too high to emit the O VII lines (0.561, 
0.568, 0.574 keV). See, e.g., figure 3 of Yoshikawa et al. 
(2003) for the oxygen line emissivity. 

It is therefore important to estimate the Galactic com- 
ponent precisely, which is possible using the offset obser- 
vation of A 1060 with Suzaku. In order to determine the 
surface brightness and the spectral shape of the Galactic 
component, we performed the simultaneous fit of the 13- 
17' and 17-27' annuli. The Galactic component is promi- 
nent in these annuli as shown in figure 6(f) and (g), how- 
ever the ICM component is still dominant almost all the 
energy range except for the OVII line. We made the si- 
multaneous fit in the whole 0.4-7.1 keV range (except 
1.825-1.840 keV), assuming one or two apec models for 
the Galactic component, and the fit results are presented 
in table 5 and figure 7. The resultant normalization of the 
apec model in table 5 is scaled so that it gives the surface 
brightness in the unit solid angle of arcmin 2 . 

It may appear that the difference between these two 
models are not large in figure 7, however, it is notable that 
the derived OVIII line intensity for the Galactic compo- 
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Fig. 7. (a) A magnification of 0.49-1.4 keV range of the BI 
(red) and FI (black) spectra in the 17-27' annulus. They are 
fitted with the apec 1 + phabs X vapec model, and the best- 
fit model are drawn by green and yellow lines. The estimated 
CXB and NXB components were subtracted, and a simultane- 
ously fit with the 13-17' annulus was conducted in the entire 
energy range of 0.4-7.1 keV. The apec 1 or vapec component 
for BI is indicated by a cyan or magenta line, respectively. 
Lower two panels show the fit residuals for BI and FI in unit 
of a. The orange dotted lines in figure 10 correspond to this 
modeling of the Galactic component, (b) Same as (a) but 
fitted with the apec 1 + apec 2 + phabs X vapec model, and the 
apec 2 component for BI is indicated by an orange line. We 
adopt this model in the spectral fit. See subsection 4.3 for 
details. 

0.17 r 1 1 , , . , , 1 . 
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0.2 0.4 0.6 0.8 

Oxygen abundance of ICM (solar) 

Fig. 8. A plot of confidence contour between kT\ (temper- 
ature of the cooler part of the two apec components) and 
the O abundance of vapec for the 17-27' annulus, in the 
simultaneous fitting of 13—17' and 17-27' annuli with the 
apec 1 + apec 2 + phabs X vapec model. The cross denotes the 
best-fit location, and the two contours represent la and 90% 
confidence ranges, from inner to outer, respectively. 

nent is by about twice larger for (a). We examined the 
improvement of the \ 2 (^X 2 = 46) with the .F-test, and 
adding the apec 2 component was justified with a large sig- 
nificance (false probability ~ 10~ 9 ). We further allowed 
the normalization of the apec 2 to become free between the 
two annuli. The derived surface brightness ratio of the 13- 
17' annulus to 17-27' was 0.9 ±0.2, hence it was consistent 
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Table 5. The best-fit parameters of the apec component(s) for the simultaneous fit of the spectra in 13-17' and 17—27' annuli. 



Fit model 


Normi * 


fcTi 


Norm 2 * kT 2 


X 7dof 






(keV) 


(keV) 


(a) apec 1 + phabs x vapec 

(b) apec 1 + apec 2 + phabs x vapec 


1.20 ±0.07 
. ... 1.66 ±0.35 


179 +u.uub 

143+ 0014 
U.14O_ 015 


0.21 ±0.05 0.737±g;g|| 


1128/994 
1082/992 



* Normalization of the apec component divided by the solid angle, H D = 7r X (20') 2 , assumed in the uniform-sky ARF calculation 
(20' radius from the optical axis of each XIS sensor), Norm = j n c n H dV / (4tt (1 + z) 2 D\ ) /Q u xlO 20 cm 5 arcmin 2 , where Da 
is the angular distance to the source. 



to be constant. We also confirmed that the 9-13' annulus 
gave the same result. We therefore concluded that the two 
apec models arc required to account for the Galactic com- 
ponent. Hereafter, the apec 1 + apec 2 +phabs x vapec model 
are utilized for the spectral fitting, otherwise stated. 

To demonstrate how sensitive the O abundance of the 
ICM is to the assumed Galactic component model, we 
present a confidence contour between fcTi (keV) of the 
apec 1 component and the O abundance (solar) of vapec 
for the outermost annulus (17-27') in figure 8. There 
appears to exist a negative correlation between the two 
parameters, because higher temperature of the Galactic 
component produces more O VIII emission line relative to 

VII, which contribute to reduce the O VIII line from the 
ICM (vapec component). Influences on the derived tem- 
perature and abundance by the modeling of the Galactic 
component will be tested in subsection 4.4, too. 

In order to take into account both existence of the 
Galactic component itself and propagation of its statistical 
error, we simultaneously fitted each annulus with the out- 
ermost annulus of 17-27'. As mentioned in the previous 
subsection, the normalization of the apec 1 + apec 2 compo- 
nent are constrained to give the same surface brightness 
between the two annuli. The temperatures of the two apec 
models were also common between the two, however their 
values (two normalizations and two temperatures) them- 
selves are left free. It is confirmed that the derived normal- 
izations, Norm i and Norm 2 , and temperatures, fcTi and 
fcT2, are consistent with the values in table 5(b) within 
the quoted errors. 

4-4- Radial Temperature & Abundance Profiles 

Before entering the spectral fit at each annulus, we in- 
vestigate the fit result in the central 0-6' region to test 
the capability in the abundance determination with the 
Suzaku XIS. This region exhibits nearly constant tem- 
perature and metal abundances (figure 10), and the con- 
tribution of the Galactic component is almost negligible 
(figure 6). The fit result is presented in table 6, in which 
most of clement abundances are allowed to be free in the 
spectral fit, except for He, C, and N, which are fixed to 

1 solar with the assumed abundance ratio of angr. 
Although the fit was not acceptable due mainly to the 

very high photon statistics than the systematic errors in 
the instrumental response, this result was useful to assess 
whether each element abundance was reasonably (sub- 
solar to ~ 2 solar) determined or not. The Ar and Ca 



abundances were reasonably determined although the er- 
rors are larger than other elements. We therefore decided 
to link the S, Ar, and Ca abundances to be the same, be- 
cause they showed similar values. The Al abundance be- 
came 0.0, which is physically strange, therefore we fixed 
the Al abundance to 1 solar. The x 2 /dof was increased to 
1898/988 with this treatment, however we confirmed that 
other parameters did not change beyond the quoted error 
range in table 6. Somehow, the Ne and Ni abundances 
were larger than other elements. This might be due to 
these clement lines could not be resolved from the Fe-L 
line complex. Note that the Ni abundance was also de- 
termined by Ni-L lines because we ignored energy range 
above 7.1 keV, while the Fe abundance was determined by 
both Fe-L and Fe-K lines. Anyway, we left these Ne and 
Ni abundances to move freely during the spectral fit. We 
also present results when the O, Ne, and Mg abundances 
are linked to have the same value in Appendix 2. 

We also found that the assumed abundance table, which 
defines the "1 solar" of each element relative to H, sig- 
nificantly affects the determination of temperature and 
abundance. In fact, difference in the abundance table 
affects the fit results in three ways. Firstly, we assume 
the Galactic component to have the apec 1 + apec 2 model 
with "1 solar". Secondly, the Galactic absorption model 
of phabs is also changed by the assumed abundance table. 
Lastly, the derived abundances are given in unit of "so- 
lar" by the vapec model. The first two effects are rather 
complicated, and details are investigated in Appendix 2. 
In this section, we treat only the solar abundance ratio 
of angr (Anders & Grevesse 1989) with the phabs absorp- 
tion model, simply because they are the default standard 
of XSPEC. In terms of the x 2 , it appears to give the min- 
imum x 2 with the lodd (Lodders 2003) abundance table 
in combination with the wabs (Morrison & McCammon 
1983) absorption model which utilizes the abundance ta- 
ble of aneb (Anders & Ebihara 1982) built-in the code 
(table 11). 

We then examined the influence of uncertainty in 
the OBF contaminant using the spectrum within 6'. 
Figure 9(a) shows a magnification of the best-fit spectrum 
in the energy range of 0.4-1.0 keV with the "nominal" 
ARF response, in which transmission of the OBF con- 
taminant was estimated according to the calibration files 
ae_xiALcontami_20060525.f its (7V = 0,1,2,3), as shown 
in figures 2 and 5(b). As described in subsection 4.2, the 
spectrum was simultaneously fitted with the 17-27' annu- 
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Table 6. Result of the spectral fit in the central region of 0—6', with the apec 1 + apec 2 +phabs X vapec model. Only the parameters 
for the vapec component are presented. The abundance of each element (except for He, C, and N, which are fixed to 1 solar with 
the assumed abundance ratio of angr) is allowed to be free in the spectral fit. Errors arc 90% confidence range of statistical errors, 
and do not include systematic errors. 



Region 


kT 
(keV) 




(solar) 


Ne 
(solar) 


Mg 
(solar) 


Al 
(solar) 


Si 
(solar) 




0-6' 


o r 4 + 0.03 


31+ 08 


-, 91 +o.io 
1 -^ 1 -o.io 


70+ ' 10 

u - ' u -0.09 


00+° 13 
u - uu -o.oo 




1871/986 




Norm * 


S 

(solar) 


Ar 
(solar) 


Ca 
(solar) 


Fe 
(solar) 


Ni 
(solar) 






302 ±3 


u - ' -0.07 




u - oo -0.16 


44+ - 17 


u - 4,3 -o.oi 


0.83t°;i 





* Normalization of the vapec component scaled with a factor of SOURCE_ratio_REG / AREA in table 4, 

Norm= SOURCE_RATIOJREG J n c n H dV / (4n (1 + z) 2 D 2 ) xl0~ 20 cm" 5 arcmin" 2 , where D A is the angular distance to the source. 



Table 7. Results of the spectral fits in the central region of 0-6' by changing the amount of the OBF contaminant. 



Contaminant O Ne Mg Al Si S Ar Ca Fe Ni 

(solar) (solar) (solar) (solar) (solar) (solar) (solar) (solar) (solar) (solar) 

nominal 0.31 L21 O70 (L00 057 057 053 044 043 083~ 

+20% 0.43 1.05 0.48 0.00 0.47 0.50 0.45 0.46 0.41 0.41 

free * 0.38 1.12 0.57 0.00 0.51 0.53 0.47 0.43 0.42 0.54 



Contaminant kT Norm'' kTi Normi X kT 2 Norm 2 t A C o 1/6 § X 2 /dof 

(keV) (keV) (keV) BI FI 

nominal 3.54 302 0.159 1.00 0.752 1.00 1.74 2.32 1878/986 

+20% 3.36 298 0.148 1.40 0.721 1.09 2.11 2.72 1746/986 

free* 3.50 301 0.159 0.96 0.765 0.95 2.18 2.38 1681/982 



* We fitted the amount of the OBF contaminant as a free parameter with the C/O number ratio fixed to 6. 
t Normalization of the vapec model, calculated in the same way with table 6 
t Ratio to the normalization of the "nominal" ARF. 

§ Column density of the OBF contaminant with chemical composition of COjyg for BI and FI in unit of 10 18 cm -2 . 



Table 8. Summary of the best-fit parameters of the vapec component for each annular region with the apec 1 + apec 2 + phabs X vapec 
model. Each annulus is simultaneous fitted with the the outermost 17-27' annulus. Errors are 90% confidence range of statistical 
errors, and do not include systematic errors. The solar abundance ratio of angr is assumed. These results are plotted in figure 10. 
The O abundance in the outer two regions becomes lower when the Galactic component is expressed by a single apec model, which 
is drawn by orange dotted lines in figure 10. The kT, O, and Mg columns at r < 6' arc slightly different when the OBF contaminant 
is increased by "+20%", which is drawn by black dotted lines in figure 10. The Ne abundance is probably not reliable because the 
Suzaku XIS cannot resolve the ionized Ne lines from the Fe-L line complex. 



Region Norm * kT 





Ne 


Mg 


Si 


S, Ar, Ca 


Fe 


Ni 


X 2 /dof 


(keV) 


(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


(solar) 



0-2' 


605 + 12 


3 


0^+0.04 





QQ+0.16 

"J 8 -0.15 


1. 


90+0.20 
zu -0.20 


0. 


.761° 


.19 
19 


0. 


,691° 


li 
.11 


0. 


.821° 


12 
.12 





40+0.03 

■ 4y -0.03 


1 


1 q+0.46 
• lo -0.45 


1240/992 


2-4' 


352 + 5 


3 


oq+o.03 

• Jt, -0.03 





9 8 +0.11 
• zo -0.12 


1. 


11+O.I6 
1J --0.15 


0. 


,65t°; 


.14 
.14 


0. 


.521° 


.08 
.08 





.521°; 


08 
.08 





/I9+0.02 
•^ z -0.02 





qf) +0.34 

yu -0.34 


1290/992 


4-6' 


198 + 3 


3 


/1O+0.04 
'^ z -0.04 





9Q+0-13 
■^-0.13 


1. 


94+0.18 
z4: -0.16 


0. 


■6718 


.15 
15 





.581°; 


09 
.09 


0. 


.531° 


.09 
.10 





^ 1+0.02 
•^ 1 -0.02 





07+0. 38 
•°'-0.37 


1345/992 


6-9' 


109 + 2 


3 


or+0.04 
• zo -0.04 





co+0.14 
• oo -0.14 


1. 


1 7+O.I6 
1( -0.15 


0. 


.591°; 


.15 
15 





.501° 


.08 
.08 


0. 


.551°; 


.09 
.09 





4f) +0.02 
•^ u -0.02 


1 


1 7+0.36 
- 1 ' -0.35 


1329/992 


9-13' 


44 + 2 


2 


S7+0.10 





oo -0.29 


0. 


oq+0.35 
oy -0.33 


0. 


■75t° 


33 
.32 


0. 


.271°; 


.17 
.17 


0. 


.191° 


20 
19 





4f) +0.05 
• w -0.05 


1 


io+0.77 
• lo -0.72 


1057/992 


13-17' 


26 + 1 


2 


r n +0.09 





79+0.29 
'' z -0.38 


0, 


44 +0.29 
**-0.28 


0. 


■78±g; 


29 
.27 





.301° 


.15 
.14 





,441° 


18 
17 





9C+0.05 
■ zd -0.05 





1R +0.60 
lo -0.16 


1082/992 


17-27' 


13 + 1 


2 


99 +0.11 
• zz -0.12 





,,0+0.41 
• 4S -0.37 





on+0.29 
JU -0.26 





.211° 


22 
.21 





,281°- 


12 
.12 





.211° 


.14 
.14 





01 +0.04 
• zl -0.04 





nf) +0.41 
UU -0.00 


_t 



* Normalization of the vapec model, calculated in the same way with table 6. 

t The 17-27' annulus was fitted simultaneously with other annuli, and the best-fit values with the 13-17' annulus are presented here. 
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Fig. 9. (a) A magnification of the best-fit spectrum at 0-6' annulus in the energy range of 0.4—1.0 keV with the 
"nominal" ARF response, in which transmission of the OBF contaminant was estimated according to the calibration files 
ae_xi./V_contami_20060525 . f its (JV = 0,1,2,3). Table 6 and the "nominal" row of table 7 show the best-fit parameters, (b) 
The best-fit spectrum with the "+20%" ARF, in which amount of the OBF contaminant was increased by +20% than the "nominal" 
ARF. The "+20%" row of table 7 shows the best-fit parameters, (c) The best-fit spectrum with the no-contaminant ARF but fitted 
with varabs X (apec 1 + apec 2 + phabs X vapec) model to consider the transmission of the OBF contaminant in the XSPEC varabs 
model with chemical composition of COyg. The "free" row of table 7 shows the best-fit parameters, (d), (e), (f) Same as (a), (b), 
(c) but for the 17-27' annulus, respectively. Each spectrum was simultaneously fitted with (a), (b), (c). 



lus shown in figure 9(d). In this figure, it is suggested 
that the absorption in the low energy band is slightly 
inconsistent between BI and FI sensors, namely, the BI 
spectrum appears to need more absorption in 0-6' where 
the amount of contamination is larger than 17-27' (see 
figure 2). We also found that the contamination mea- 
surement by RXJ1856. 5-3754 on October 25, 2005 in the 
XIS hardware paper (figure 14 of Koyama et al. 2006) 
indicates by about 20% larger amount of the OBF con- 
taminant than the "nominal" value for BI (XIS1). 

We therefore generated ARF responses changing the 
amount of the OBF contaminant by +20%. The fit re- 
sult is shown in figure 9(b), (e) and table 7 ("+20%" 
row). Note that the CXB background to subtract was 
also modified in this fit. The fit residual for 0-6' was 
improved, and the x 2 was decreased by A^ 2 = 132. We 
further tested the fit by adding an absorption of the OBF 
contaminant with chemical composition of COi/g to the 
fit model using the XSPEC varabs model. 2 The fit re- 
sult is shown in figure 9(c) and (f). The best-fit values 
are summarized in table 7 ("free" row), and the \ 2 was 
improved by A% 2 = 197. The derived amount of contam- 
inant for BI was by 25% larger than the "nominal" value, 
and by 3% larger for FI. We think that the deviation 

2 In this method, we used the same CXB background with the 
"nominal" ARF, and the surface brightness of the Galactic com- 
ponent was slightly different between 0—6' and 17—27'. 



from the "nominal" value for BI is beyond the calibration 
uncertainty, however, almost all of the best-fit values at 
the "free" row in table 7 are between values at "nominal" 
and "+20%". We therefore consider the +20% result as 
the systematic error range due to the uncertainty in the 
OBF contaminant. 

Considering these systematics, results of the spectral fit 
at each annulus are summarized in table 8 and figure 10. 
We tested the results by changing the background nor- 
malization by ±10%, and they are plotted in green dotted 
lines in figure 10. The systematic error due to the back- 
ground estimation is almost negligible. Difference in the 
best-fit values by modeling the Galactic component with 
a single apec was investigated, and they were indicated by 
orange dotted lines. The differences are within the sta- 
tistical error for the inner five annuli, however, O abun- 
dance becomes lower than the 90% confidence error at 
the outer two annuli as indicated in figure 7, and temper- 
ature and Fe abundance become lower at the outermost 
annulus. This is due mainly to the fact that the XIS can- 
not resolve the ICM O VIII line from the Galactic O VIII 
by redshift. The systematic error range due to the un- 
certainty in the OBF contaminant is indicated by black 
dotted lines. It is sometimes larger than the statistical er- 
rors at small radii (r < 6'), particularly for kT ', O, and Mg 
abundances. Though Ne abundance is significantly larger 
than other elements at small radii (r < 9'), it is probably 
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Fig. 10. (a) Radial temperature profiles derived from spectral fit of the Suzaku (black) and XMM-Newton (red) spectra at each 
annulus. The horizontal axis denotes the projected radius and deprojection are not conducted. The same data with table 8 are 
used for Suzaku. The black dotted lines correspond to shifts of the best-fit values by changing thickness of the OBF contaminant 
by +20%. The green dotted lines denote those when the estimated CXB and NXB levels are changed by ±10%. The orange dotted 
line shows the best-fit value when the Galactic component is modeled by a single apec. Regarding XMM-Newton, the phabs X vapec 
model are used for the spectral fit, and the O, Ne, Mg, and Ni abundances are fixed to 1 solar in the assumed abundance table of 
angr. (b)— (h) Radial abundance profiles derived and plotted in the same way as (a). 



not reliable because the Suzaku XIS cannot resolve the 
ionized Ne lines from the Fe-L line complex. 

We also plot the XMM-Newton results for the inner 
five annuli with red diamonds, in which the phabs x vapec 
model are used for the spectral fit, and the Ne, Ar, Ca and 
Ni abundances are fixed to 1 solar. Namely, the Galactic 
component is ignored. They are almost consistent with 
the Suzaku results, although XMM-Newton appears to 
give slightly higher abundance for Fe and Si, while lower 
abundance for S. 

4--5. Direct Comparison of O VII and O VIII Intensities 

We also investigated the surface brightness of the O VII 
and O VIII emission lines in order to estimate the O abun- 
dance of the ICM directly from the line intensities. The 



surface brightness of the O VII and O VIII were derived by 
fitting the annular spectrum with a power-law + gaussian 
+ gaussian model. In the fitting, we fixed the Gaussian a 
to be 0, and allowed the energy center of the two Gaussians 
to vary within 555-573 eV or 648-658 eV for OVII or 
O VIII, respectively. The derived line intensities are sum- 
marized in table 9 and figure 11(a). There is a clear excess 
of the OVIII intensity towards the cluster center, while 
OVII is consistent with constant. This is a strong evi- 
dence that the O VIII emission is associated with the ICM 
itself, on the other hand, OVII might be due mainly to 
the Galactic origin. 

We compared the surface brightness of the oxygen lines 
with that at a neighbor of A 1060, the NGC 2992 field. 
Though NGC 2992 is located at 18° offset from A 1060, 
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Fig. 11. (a) Line intensities of Ovu and Oviii at each annulus of A 1060 in unit of photons cm -2 s" 1 sr -1 . Each value is shown 
in table 9. (b) The O VIII emission line intensities in table 9 divided by the surface brightness of the 2—5 keV continuum are plotted 
in green diamonds against the radius of each annulus. The red or black diamonds correspond to the ratios when the O VIII intensities 
of the NGC 2992 field or the outermost annulus (17-27') are subtracted as an offset. Dashed lines represent expected ratios for the 
vapec model with 0.25, 0.5, and 1 solar abundances, assuming the observed temperature at each annulus (table 8). 



Table 9. Line intensities of Ovu and Oviii at each an- 
nulus of A 1060 and the NGC 2992 field in unit of pho- 
tons cm -2 s — 1 sr — x . These intensities are derived from the 
spectral fit with power-law + gaussian + gaussian model (see 
figure 12 for the NGC 2992 field), assuming the uniform-sky 
ARF response. Intensities by McCammon et al. (2002) mea- 
sured with a high resolution microcalorimeter array for a large 
sky area of ~ 1 sr are also presented. 



Region O VII O VIII 

A 1060 center (0-2') s> 17.1±9.2 

A 1060 center (2-4') 7.2±7.2 12.1±5.0 

A 1060 center (4-6') 10.0±4.4 8.6±3.3 

A 1060 center (6-9') 7.3±2.0 6.9±1.4 

A 1060 offset (9-13') 8.3±2.8 5.7±1.8 

A 1060 offset (13-17') 8.6±1.7 3.5±0.9 

A 1060 offset (17-27') 8.6±1.2 2.0±0.5 

NGC 2992 field 4.1±0.6 0.6±0.2 

McCammon et al. (2002) 4.8±0.8 1.6±0.4 



o.i 



the Galactic latitude of A 1060 and NGC 2992 is sim- 
ilar (tablel), and the neutral hydrogen column density 
(A r H = 5.1 x 10 20 cm~ 2 ; Dickey & Lockman 1990) is com- 
parable to A 1060. NGC 2992 [z = 0.007710) has an active 
galactic nuclei, which is expected to emit no oxygen lines. 
In addition, we excluded the main target with a radius of 
5' centered on NGC 2992, and utilized the outer region 
for the spectral analysis in figure 12. 

The surface brightness of the O VII and O VIII emissions 
in the outermost annulus of the A 1060 observation was 
roughly twice larger than the NGC 2992 field. Line inten- 
sities by McCammon et al. (2002) measured with a high 
resolution microcalorimeter array for a large sky area of 
~ 1 sr are also presented in table 9, and the O VIII inten- 
sity was consistent with the outermost annulus (17-27') 
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Fig. 12. The 0.46-0.74 keV spectra of the NGC 2992 field 
with BI (red) and FI (black) are fitted with the power-law + 
gaussian + gaussian model. The red and black lines represent 
the best-fit model, and the model components for BI are plot- 
ted by green (power-law) and blue (two gaussian) lines. The 
central energies of two Gaussians are constrained to be within 
±5 eV range of O VII and O VIII lines, and <r of Gaussians are 
fixed to 0.0. 

of the A 1060. It is suggested that both the O VII and and 
O VIII lines in our observation are somewhat mixture of 
the ICM and the Galactic origin. However, we could not 
confirm excess O VII emission towards the cluster center, 
we concentrate on the O VIII emission. 

Then we calculated the surface brightness ratio of the 
O VIII line divided by the 2-5 keV continuum at each an- 
nulus, as presented in figure 11(b) by green diamonds. 
Here, we subtracted the O VIII intensity at the NGC 2992 
field (red diamonds) or the intensity at the outermost an- 
nulus (black diamonds) , as an "offset" due to the Galactic 
component. The dashed lines represent the calculated ra- 
tio by the XSPEC simulation fixing the O abundance of 
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ICM to 0.25, 0.5, and 1 solar. The ICM temperature at 
each annulus is assumed to have the measured value in fig- 
ure 10(a). Actual value of the O abundance is supposed 
to lie between the black and green diamonds. The radial 
O abundance profile obtained in figure 10(b) is mostly 
consistent with this plot. 

4-6. Central Cool Component of A 1060 

Chandra observation resolved two central elliptical 
galaxies of A 1060, NGC 3311 and NGC 3309 (Yamasaki 
et al. 2002). We estimated the flux of an additional cool 
component for the central galaxies with the vapec model 
of kT = 0.8 keV and 0.5 solar abundance at z = 0.0114, 
and obtained an upper limit (90% confidence level) to be 
L x = 5.2 x 10 40 erg s -1 (0.4-4.5 keV) in the central region 
within r < 2'. In this fit, temperatures of the Galactic 
components (two apec) were fixed at 0.146 and 0.662 keV, 
and only the ICM component was allowed to be free. This 
upper limit is consistent with the sum of the ISM flux 
of the two elliptical galaxies by Yamasaki et al. (2002), 
Lx = (1-8 ± 0.6) x 10 40 erg s' 1 (0.4-4.5 keV). 

5. Discussion 

5.1. Temperature Profile 

The temperature in A1060 observed with Suzaku shows 
a decline from ~ 3.4 keV in the central region (r ^ 5' ~ 
70 kpc) to ~ 2.2 keV in the outermost annulus of 17-27' 
- 240-380 kpc - 0.16-0.25 n 80 . This feature is roughly 
consistent with the previous result for the inner region (r < 
13') measured with XMM-Ncwton (Hayakawa et al. 2006), 
although our Suzaku result indicate slight flattening of the 
temperature profile in the outer annuli. Hayakawa et al. 
(2006) suggested a faster decline of the temperature pro- 
file of A 1060 than the average T/ (T) curve by Markevitch 
et al. (1998) obtained from 30 nearby clusters (exclud- 
ing A 1060) with ASCA observations. However, now it 
appears to become consistent with them by taking the 
emission-weighted average temperature, k(T) — 2.5 keV. 

Recent measurements of 13 nearby relaxed clusters (ex- 
cluding A 1060) with Chandra by Vikhlinin et al. (2005) 
showed that the temperature reaches at a peak at r ~ 
0.15 riso and then declines to ~ 0.5 of its peak value at 
r ~ 0.5 riso) in good agreement with Markevitch ct al. 
(1998). They also found that clusters whose tempera- 
ture profiles peak at r < 70 kpc < 0.15 ri 80 (MKW4 and 
RX J1159+5531) shows a larger peak-temperature-ratio of 
T p /(T) - 1.35. This is exactly the case of A 1060, which 
shows the peak-temperature of kT p = 3.42 ± 0.04 keV at 
r = 4-6' ~ 70 kpc with our Suzaku measurement. Using 
the relation of T p /(T) <~ 1.35, the emission-weighted av- 
erage temperature is calculated to be k(T) = 2.5 keV, as 
suggested in the previous paragraph. These three corre- 
spond to clusters without cooling core, hence this appar- 
ent steep temperature decline at r < 0.2 ri 80 might be a 
common feature for the non-cooling flow clusters. 

In order to determine actual value of k{T) for A 1060, 
further offset observations toward ri 80 ~ 110' are required. 
Since the wide and energy-dependent point-spread func- 



tion of ASCA tended to give a flatter temperature profile 
for A 1060 (Tamura et al. 1996; Furusho ct al. 2001), stud- 
ies of the T I (T) curve in cluster outskirts with Suzaku will 
be useful to look at the true temperature features. 

Dolag et al. (2004) conducted a smoothed particle hy- 
drodynamics (SPH) simulation to reconstruct the struc- 
ture of the Local Universe, which well reproduces the posi- 
tion, mass, and temperature of A 1060 ("Hydra" of table 1 
in Yoshikawa et al. 2004). The simulated temperature of 
A 1060 was 3.4 keV within r < 20' - 0.18 n 80 . The radial 
temperature profile of simulated A 1060 is presented by 
Kawahara et al. (2007), based on the updated version of 
the SPH simulation by Dolag et al. (2005). It reproduces 
the shape of the temperature profile obtained with our 
Suzaku observation pretty well, although the temperature 
in the simulation is higher than our observation. This is 
partly because our temperature profile in figure 10(a) is 
plotted against the projected radius, in which the ICM 
emission is integrated along the line of sight, while the 
simulated temperature is plotted against the spherical ra- 
dius. It is pointed out that the spectroscopic temperature 
in the real X-ray observation tend to give lower temper- 
ature than the emission-weighted temperature calculated 
in the simulation due to the multi-phase nature of the 
ICM not only from the radial profiles but also from the 
local inhomogeneities (Mazzotta et al. 2004; Kawahara et 
al. 2007). 

There is a small flat-top in the observed temperature 
profile within r < 5' ~ 70 kpc in figure 10(a). This feature 
might indicate an initial phase in the formation of cooling 
core. The SPH simulation predicts that the period of the 
last major merger occurred at an age of ~ 3 Gyr ago (z ~ 
0.3). This implies that the temperature flat-top has been 
made within recent ~ 3 Gyr. The thermal conduction 
length scale is calculated as 



?*cond — 155 



^cond 



3 Gyr) 



0.01 cm- 



kT 
3kcV 



kpc, 



where i CO nd is the thermal conduction time, n c is the elec- 
tron density, assuming the classical Spitzer conductivity. 
This length scale is similar to the size of the flat-top. 

In subsection 4.6, only an upper limit of 5 x 10 40 erg s _1 
for the cool component in the ICM was obtained. It con- 
firms that the sharp cusp structure in the Chandra image 
of A1060 (Hayakawa et al. 2004) really reflects the gravi- 
tational mass concentration and does not come from cool 
gas of high emissivity. 

5.2. Abundance Profiles 

In this paper, we obtained emissivity weighted abun- 
dance profiles of Ni, Fe, Si, S, Mg, Ne and O in the ICM 
of A 1060 up to a radius of 27' ~ 380 kpc, although Ni 
and Ne abundances are not reliable due to the strong and 
complex Fe-L line emissions. Abundances of Si, S, and 
Fe decrease from ~ 0.7, 0.8 and 0.5 solar in the central 
region, to ~ 0.3, 0.2 and 0.2 solar, respectively, in the out- 
skirts of the cluster. De Grandi & Molcndi (2001) showed 
that non-cooling flow clusters do not exhibit a steep abun- 
dance gradient due possibly to disruption of the central 
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Fig. 13. (a) Abundance ratios of O (black), Mg (red) and Si (blue) divided by Fe. (b) Abundance ratios of O, Ne, Mg, Si, and S, 
divided by Fe for Fornax cluster, NGC 1404 (Matsushita et al. 2007a), NGC 720 (Tawara et al. 2007), and A 1060. In both plots, 
the solar abundance ratio of angr (Anders & Grevesse 1989) is assumed, and only the statistical errors are plotted. 
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Fig. 14. (a) Cumulative Fe (red) and O (black) mass profiles of A 1060. (b) Integrated iron mass-to-light ratio (IMLR) and oxygen 
mass-to-light ratio (OMLR) in unit of Mq/Lq in the blue band, (c) Differential IMLR and OMLR calculated at each annulus. The 
annular range of 2—9' are not plotted because there are no member galaxies cataloged by Christlein & Zabludoff (2003). See text 
subsection 5.3 for details. 



cool cores after major merger events. A 1060 shows clear 
gradients for Fe, Si, and S, although it is a very relaxed 
system without the cD galaxy and cooling core. 

The other elements, Mg and O, show somewhat flatter 
radial distributions. The measured Mg abundance is con- 
sistent with constant up to the 13-17' annulus, namely 
out to 240 kpc from the cluster center. Oxygen also ex- 
hibits a flat distribution or a slight increase at r ~ 10'. We 
have to note that the O abundances depend on the esti- 
mation of the Galactic component as described in subsec- 
tion 4.3. The O abundance in the outer annuli is partic- 
ularly affected by the surface brightness and temperature 
of the cooler Galactic component (kT% ~ 0.15 keV). The 
O and Mg abundances at r < 6' are also affected by the 
uncertainties in the OBF contaminant (subsection 4.4). 
We have tried to take into account all the possible sys- 
tematic errors, and the results are presented in figure 10. 
In addition, systematic effects caused by the definition of 
the solar abundances were addressed and described in ap- 
pendix 2. 

In order to compare the relative variation in the abun- 



dance profiles, we show abundance ratios of O, Mg, and 
Si divided by Fe as a function of projected radius in fig- 
ure 13(a). Apparently, the profiles of Si/Fe, O/Fe, and 
Mg/Fe show different gradients. The Si/Fe ratio is con- 
sistent with a constant value, while the O/Fe and Mg/Fe 
ratios seem to increase with radius excluding the outer- 
most annulus (17-27'). This trend is similar when the O, 
Ne, and Mg abundances are constrained to have the same 
values in different solar units in figure 19(c). 

Recent Suzaku observations also showed the O and 
Mg abundances in poorer systems: an elliptical galaxy 
NGC 720 (Tawara et al. 2007), the Fornax cluster, and 
NGC 1404 (Matsushita et al. 2007a). The Fe ratios of ele- 
ments for these systems are compared in figure 13(b). The 
ratios Si/Fe and S/Fe are almost the same with A 1060 
in other systems. Somehow, the Mg/Fe ratio in A 1060 
is ~ 50% higher than those in other systems. However, if 
one considers the influence of the OBF contaminant, this 
ratio can drop and the difference becomes insignificant. 
Abundance ratio of O/Fe shows almost the same values 
with A 1060 at their centers, and in the outer regions 



No. ] Temperature and Abundance Profiles of Abell 1060 with Suzaku 15 

Table 10. Comparison of IMLR and OMLR with other systems. 



IMLR OMLR < r kT Reference 

NGC 720 1 x 10~ 4 4 x 10~ 4 25 kpc - 0.56 keV Tawara et al. (2007) 

Fornax 4 x 10~ 4 2 x 10~ 3 130 kpc - 1.3 keV Matsushita et al. (2007a) 

Ccntaurus 4 x 10~ 3 3 x 10~ 2 190 kpc ~ 4 keV Matsushita et al. (2007b) 

A 1060 4.0 x lO" 3 4.3 x 10~ 2 380 kpc - 3 keV This work 



A 1060 indicates a higher value than others. 

This peculiar behavior of the O/Fe ratio has been re- 
ported by Tamura et al. (2004) with XMM-Newton as an 
average of 19 clusters which are mainly composed of X-ray 
bright and relaxed clusters with a cD galaxy (excluding 
A 1060). The O/Fc ratio is ~ 0.7 solar within r < 50 kpc, 
which increases to ~ 1.5 solar at r ~ 100 kpc, and more at 
r > 200 kpc, in figure 5(b) of Tamura et al. (2004), assum- 
ing the solar abundance ratio of angr (Anders & Grevesse 
1989). This is much alike the O/Fe plot in figure 13(a). 

Since these differences are not simply explained by the 
ion mass, we need to consider that the enrichment pro- 
cesses have a significant difference between O and Fe, for 
example. It is generally considered that enrichment of 
O and other SN II originated metals has occurred in the 
early stage of cluster formation, certainly before the last 
merger epoch. The relatively flat distributions of O and 
Mg indicate that the early metal enrichment has caused 
these features. One plausible explanation is the enrich- 
ment in the form of starburst-driven galactic winds (e.g. 
Strickland & Stevens 2000). Recent numerical simula- 
tion of the cluster metal enrichment indicates that even 
a metallicity peak in an intermediate radius can be cre- 
ated by enrichment through galactic winds (Kcpferer et al. 
2006, priv. comm.), which resembles our oxygen feature. 

The other possibility is the early metal enrichment by 
massive Population III stars (e.g. Matteucci et al. 2006) 
In this case, a large fraction of the intergalactic space 
would be enriched with metals. One big riddle in the ob- 
tained abundance profiles is that the Mg abundance at the 
central region appears to be significantly higher than the 
O abundance, although both O and Mg are mainly pro- 
duced by SN II. Loewenstein (2001) suggests that preen- 
richment by a generation of massive Population III stars 
may account for the abundance anomalies of O, Si, and 
Fe observed in the ICM. However, Yoshida et al. (2004) 
concluded that metals originating from the earliest gener- 
ation of stars cannot be responsible for the observed abun- 
dance anomalies in the ICM. Future quantitative studies 
of metal abundances in the outer regions (~ riso) or even 
outside of clusters will be important to locate the origin of 
metals. Instruments with much higher energy resolution 
such as microcalorimeters will be desired for these studies. 

5.3. IMLR and OMLR 

Combining the X-ray luminous gas mass profile by 
Hayakawa et al. (2006) with XMM-Newton and the abun- 
dance profiles obtained with Suzaku, we calculated the 
cumulative iron or oxygen mass profiles in figure 14(a). 



Strictly speaking, the abundance profile is derived from 
the projected spectrum along the line of sight, we ap- 
proximated it to the spherical distribution. The iron 
and oxygen mass within r < 380 kpc is calculated to be 
~ 1.4 x 10 9 M G and 1.6 x 10 10 M Q , respectively. 

Based on the member galaxy catalog (69 galaxies within 
the projected radius, r < 27' ~ 380 kpc) by Christlein & 
Zabludoff (2003), we also calculated the iron mass-to- light 
ratio (IMLR) and oxygen mass-to-light ratio (OMLR) 
in figure 14(b) and (c). We utilized the observed red- 
shift to estimated the 3-dimensional distribution of the 
galaxy, assuming the spherical symmetry. The IMLR and 
OMLR are calculated to be ~ 4.0 x 10~ 3 M Q /L Q and 
4.3 x 10~ 2 Mq/Lq within r < 380 kpc. Both are sub- 
jected to large errors, and we tentatively adopted factors 
of two and three, respectively, in figure 14. 

We summarize recent measurements of IMLR and 
OMLR in table 10. NGC 720 and the Centaurus cluster 
exhibit similar values with A 1060 when compared at the 
same radius, while the Fornax cluster shows significantly 
smaller values by an order of magnitude. It is suggested 
that the smaller systems with lower gas temperature tend 
to show lower IMLR in Makishima et al. (2001), however 
its deviation is outstanding. 

Both the integrated IMLR and OMLR steeply increase 
up to ^100 kpc and seem to reach almost maximum at 
100-200 kpc. This trend is similar to the IMLR obtained 
with XMM-Newton for M 87 and the Centaurus cluster by 
Matsushita et al. (2007b) within r < 100 kpc. The steep 
increase in the r < 100 kpc region suggests that the Fe 
and O ions which was synthesized in the central galaxies 
have diffused to the ICM. On the other hand, as pointed 
out in Ezawa et al. (1997), abundance gradient over a few 
hundred kpc scale in clusters of galaxies should follow the 
mass ratio between the galaxies and the ICM gas. This is 
because heavy ions released from galaxies without much 
kinetic energy would not diffuse out more than 10 kpc 
over the Hubble time. (Bohringer et al. 2004) indicated 
that it takes about 10 10 yr to synthesize Fe mass within 
r < 100 kpc range. 

These considerations suggest that both Fe and O en- 
riched gas was released to the ICM space with significant 
kinetic energy and/or the gas stripping was efficiently oc- 
curred in the ^100 kpc range due to the galactic motion. 
We also note that two central galaxies (NGC 3311 and 
NGC 3309) of A 1060 are slightly off-center (~ 0.3' and 
2') from the X-ray peak of the cluster (see figure 1(b) of 
Hayakawa et al. 2006), and the line of sight velocity is cal- 
culated to be 180 km s -1 and 660 km s^ 1 . Ram-pressure 
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stripping is also likely to be the reason why these central 
galaxies show very compact X-ray halo (Yamasaki et al. 
2002). 

The obtained O mass in figure 14(a) provides the rough 
estimate of the SN II rate. Theoretical calculations 
(Nomoto et al. 2006; Thielemann et al. 1996; Tsujimoto 
et al. 1995) predict the O product of - 1.5 M Q for 20 M Q 
progenitor mass and ~ 9 M Q for 40 M Q progenitor. This 
requires at least ~3x 10 9 SN II in the whole cluster since 
the formation. 

5.4- Warm Component in the Galactic Emission 

As shown in figure 11(a), the association of O VIII line 
and the ICM is evident. The surface brightness of OVII 
lines are brighter than the level of NGC 2992 region or 
than the average of wider sky region (McCammon et al. 
2002), but consistent with the uniform distribution over 
the cluster. A large fraction of O VII photons is considered 
to be originated from the Milky Way Halo, whose typical 
temperature measured with XMM-Newton is kT ~ 0.2 
keV (Lumb et al. 2002). If there is a cooler component 
like Local Hot Bubble, kT ~ 0.08 keV by Lumb ct al. 
(2002); Snowdcn ct al. (1998), or Warm-Hot Intergalactic 
Medium, it is hard to detect it under the strong emission 
from the ICM. Our spectral fit showed that a fcT-0.7 keV 
component was required to reproduce the spectral struc- 
ture around the Fe L-line complex. The surface bright- 
ness of this component shows a flat distribution in our 
field of view (subsection 4.3). The ROSAT 3/4 keV image 
shows an enhancement with an angular extent of ~ 10 x 20 
deg 2 , which may be an extension of the North Polar Spur 
(Snowdcn ct al. 1997). These features suggest that the 
^0.7 keV component is a Galactic origin, but a further 
study is certainly needed and it is beyond the scope of 
this paper. 
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Appendix 1. Fraction from Corresponding Sky 

Since angular resolution of X-Ray Telescopes (XRT; 
Scrlemitsos et al. 2007) of Suzaku (~ 2' in half power di- 
ameter; HPD) is not as good as those of Chandra (~ 0.5") 
and XMM-Newton (~ 15"), fraction of photons may be- 
come significant that come from outside of the sky direc- 
tion corresponding to the extraction region of spectra on 
the XIS detector plane. This was a severe problem for the 
ASCA GIS (Ohashi ct al. 1996; Makishima et al. 1996), 
which employs similar type of thin-foil-nested reflectors 
for the XRT (Serlemitsos et al. 1995) with - 3.6' HPD, 
and the fraction sometimes became more than half of to- 
tal detected photons as seen in figure 4 of Kikuchi ct al. 
(1999) for the MKW 3s and 2A 0335+096 clusters. It 



made the cluster analysis further complicated that energy 
dependence of the fraction could not be neglected with 
the ASCA XRT/GIS system. 

The situation has been much improved for Suzaku than 
ASCA. The Suzaku XRT/XIS system adopts longer focal 
length of 4.75 m than ASCA (3.5 m) and the sky cover- 
age of the XIS detector (18' x 18') is smaller than the GIS 
(~ 22' radius), so that the energy-dependent vignetting ef- 
fect is much smaller. Moreover, improvement of the mirror 
surface has reduced scatters on the reflector and equip- 
ment of the pre-collimator dramatically decreased the so- 
called "stray-lights" mainly from outside of the XIS field 
of view (Serlemitsos et al. 2007). 3 

We therefore estimated the fraction of photons out- 
side the extraction region using a simulator of the Suzaku 
XRT/XIS system, "xissim" version 2006-08-26 (Ishisaki et 
al. 2007) , in the following way: (1) Assuming the double- 
(3 model surface brightness profile (table 3), 2,000,000 
count of monochromatic incident-photons (in 0.5, 1, 2, 
4, or 8 keV) were generated by a "mkphlist" task. 
(2) Simulated event files were created for both the central 
and offset observations using the "xissim" task, with the 
xis_efficiency parameter set to "no" to save photon statis- 
tics. (3) Using the RA and DEC columns in the simulated 
event files, the event files were splitted into 7+1 sky re- 
gions corresponding to the extraction annuli of 0-2', 2-4', 
4-6', 6-9', 9-13', 13-17', and 17-27', plus outside of them, 
r>27'. 4 (4) The detected photon count at each extrac- 
tion region on the XIS detector plane was calculated for 
each splitted event file using the "xselect" task. 

In this way, we estimated the fraction of photons for the 
7 extraction annuli from the 8 sky regions in 5 energies 
for the 4 XIS sensors. We plot examples of BI (XIS1) in 
figure 15. It is confirmed that more than 65% of photons 
are coming from the corresponding sky directions by the 
"xissim" simulations. The energy dependence is almost 
negligible within ~ 5% in figure 15(h). Fractions from the 
second next to the noticed annulus are less than ~ 1% at 
the offset observation, even though the surface brightness 
of the central part of the cluster is much brighter by orders 
of magnitude (figure 4). This is primarily owing to the 
existence of the pre-collimator for the Suzaku XRT. These 
results are almost the same with the other XIS sensors. 
Note that this simulation docs not include the CXB and 
the Galactic background, which are almost uniform on the 
sky so that the spatial response is different from figure 15. 

Appendix 2. Dependence on Abundance Tables 

There are several abundance tables usable in XSPEC, 
such as angr (Anders & Grevesse 1989), aneb (Anders & 
Ebihara 1982), feld (Feldman 1992), and lodd (Lodders 
2003). The elemental abundances relative to hydrogen in 

3 The "stray-lights" reach to the focal plain through abnormal 
paths, e.g., reflection by only the secondary mirror, or the mirror 
backside reflection (Mori et al. 2005). 

4 The "fselect" task with the angsep function was utilized to calcu- 
lated the angular separation of events from the assumed center 
of the cluster, (RA, Dec) = (10 h 36 m 42F8, -27°31'42"). 
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Fig. 15. (a)-(g) Fraction of photons detected in the 0-2' to 17-27' extraction regions on the BI (XIS1) sensor plotted against the 
original sky directions of incidence, which is estimated by the "xissim" simulation. Different energies are plotted in different colors, 
(h) Photon fraction of each annulus coming from the corresponding sky region, i.e., maximum values of panels (a)-(g) are plotted. 
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Fig. 16. (a) The assumed metal abundances relative to H in number by the angr, lodd, feld, and aneb tables of XSPEC vll.3.0t. (b) 
Same as (a) but normalized by the angr abundances, (c) Comparison of absorption by the phabs model with angr abundance table, 
phabs with lodd, and the wabs model in XSPEC. The wabs model assumes the abundance ratio of aneb built-in the code. Neutral 
hydrogen column density of TVjj = 4.9 X 10 20 cm -2 and "1 solar" abundance, which is different among the tables, are assumed. The 
photoelectric absorption cross-section of beme (Balucinska-Church & McCammon 1992; Yan et al. 1998) is used. 



number denned by these tables are plotted in figure 16(a), 
and ratios of each element to angr is shown in figure 16(b). 
Although feld is not used in our spectral analysis, it is 
added to compare recent Suzaku results (Matsushita et 
al. 2007a; Tawara et al. 2007) which adopt this table. As 
seen in figure 16(b), angr and feld are essentially the same 
except for Fe abundance, which is 1.44 times larger for 
angr than feld. 

The angr table is the de facto of XSPEC, however, Fe 
abundance is significantly different from others. This is 
because angr is based on the solar photospheric measure- 
ment of Fel (Blackwell et al. 1984), which disagrees with 
the meteorites measurement adopted in aneb. Discussions 
afterwards seem to prefer the lower abundance by mete- 
orites, because Fe abundance determined by Fell, which is 
dominant in the solar photosphere, is consistent with the 
meteorites value (Biemont et al. 1991; Raassen & Uylings 
1998). 

It is notable that C, N, O, and Ne abundances in lodd 



is significantly lower (~ 60% of angr). These values are 
based on the three-dimensional time-dependent hydrody- 
namical model solar atmosphere, in which departures from 
local thermodynamical equilibrium (LTE) and blend of 
Ni I line with the forbidden [O i] A6300 are considered 
(Allende Prieto et al. 2001; Allende Prieto et al. 2002). 
The lodd table in XSPEC is the recommended elemental 
abundance of the solar photosphere (table 1 in Lodders 
2003), however, he states that the photo-solar abundances 
as representative of the solar system should be slightly 
larger (by ~ 1.2 times for C, N, and O) for helium and 
heavy-element due to the settling effects. 

The decrease of these abundant elements in the Inter 
Stellar Medium (ISM) also leads to the difference in the 
assumed Galactic absorption as seen in figure 16(c). The 
phabs absorption model with lodd abundance table gives 
by ~ 15% smaller absorption at the maximum in our no- 
ticed energy range of 0.4-7.1 keV in the spectral fitting, 
while the difference between wabs and phabs with angr is 
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Fig. 17. Same as figure 10 of the Suzaku result for the black crosses (phabs angr). Other fit results are also plotted in blue, red, 
and green diamonds, when different abundance ratio (lodd) and/or different absorption model (wabs) are assumed. 



Table 11. List of x 2 /dof for each fit. 
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negligible. This is because we have assumed the neutral 
hydrogen column density of Nr. — 4.9 x 10 20 cm" 2 (Dickey 
& Lockman 1990) with "1 solar" abundance for the ISM. 
In fact, the ISM consists of monoatomic gas, molecules, 
and grains, so that it is suggested to have smaller abun- 
dance than the solar abundance by Wilms et al. (2000). 
Shaver et al. (1983) have reported a radial abundance gra- 



dient in the Milky Way Galaxy of ~ 16% per kpc. 

Considering these uncertainties, we have tested the 
spectral fit of A 1060 with four kinds of combinations 
in the absorption model and the abundance table, phabs 
angr, wabs angr, wabs lodd, and phabs lodd. We further 
tested the fit when the abundances of O, Ne, and Mg are 
linked to have the same value relative to the assumed "so- 
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Fig. 18. (a) O/Fe ratio, (b) Mg/Fe ratio, and (c) Si/Fe ratio, which are same as figure 13(a) for the black crosses (phabs angr) 
in these figures (a)— (c). Other fit results are also plotted in blue, red, and green diamonds, when different abundance ratio (lodd) 
and/or different absorption model (wabs) are assumed. 
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Fig. 19. Radial profiles of (a) temperature, (b) O = Ne = Mg abundance, and (c) abundance ratio to Fe, when the abundances of 
O, Ne, and Mg are constrained to have the same value for the assumed abundance tables of angr or lodd. Meaning of the markers 
are same as figures 17 and 18. 



lar" abundance, in two reasons: they are all supposed to 
be primarily the SN II products; O and Ne abundances 
are susceptible to the Galactic component in outer annuli. 
Note that the O : Ne : Mg ratios are significantly different 
between angr (1 : 0.145 : 0.045) and lodd (1 : 0.151 : 0.072). 

Results are summarized in figures 17-19 and a list of 
X 2 /dof is presented in table 11. As expected from fig- 
ure 16(c), the phabs angr and wabs angr results are almost 
identical. The phabs lodd gives slightly larger tempera- 
ture than others, which leads the O abundance profile in 
figure 17(b) to be different behavior from others. Other 
discrepancies in the derived abundances can be explained 
by difference in the definition of "1 solar". 

Since y 2 values for phabs lodd are generally larger than 
other models, the Galactic absorption model by phabs lodd 
is considered to be inappropriate with our X-ray data. On 
the other hand, wabs lodd model shows the minimum \ 2 m 
table 11. This is interesting because the wabs absorption 
model assumes aneb abundance table built-in the code. 
This may indicates that our modeling of the Galactic com- 
ponent is too naive, and/or the Galactic abundance gra- 
dient is responsible so that effective elemental abundance 
between absorption and emission arc different. 

It is notable that shape of the radial abundance profiles 
are quite similar among these four models in figure 19 
when O = Ne = Mg is assumed. Difference from Fe or Si 



is evident as seen in figures 18(c) and 19(c): Fe and Si 
abundances decreases with radius with roughly constant 
Si/Fe ratio, while O = Ne = Mg are almost constant within 
r<17'. The higher (0 = Ne = Mg)/Fe ratio than "1 solar" 
is also suggested for both angr and lodd abundance tables. 
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